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Through the analysis of HELIOS-3 data obtained at the 
CERN SPS, we demonstrate the importance of secondary pro- 
cesses for dilepton production in heavy-ion collisions in the 
intermediate invariant mass region. We find that while the 
dilepton spectra between 1 to 2.5 GeV from proton-induced 
reactions can be attributed to the decay of primary vector 
mesons, charmed hadrons, and initial Drell-Yan processes; the 
strong enhancement seen in the heavy-ion data as compared 
to the background comes mainly from the secondary processes 
which are germane to heavy-ion collisions. Furthermore, we 
find nai — > 11 to be the most important process in this mass 
region, as was found by thermal rate calculations. We em- 
phasize the constraints on the elementary cross sections by 
the experimental data from e + e~ annihilation. 
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The experimental measurement and theoretical inves- 
tigation of dilepton production in nuclear collisions con- 
stitutes one of the most active and exciting fields in 
heavy- ion physics Q|. Because of their relatively weak 
final-state interactions with the hadronic environment, 
dileptons are considered ideal probes of the early stage 
of heavy-ion collisions, where quark-gluon-plasma (QGP) 
formation and chiral symmetry restoration are expected 

11; 

Dilepton mass spectra in heavy ion collisions can be 
divided into three regions. The region below (~ 1 
GeV) is dominated by hadronic interactions and hadronic 
decays at freeze-out. In the intermediate- mass region, 
rrirj, < M < to j /y, the contribution from the thermalized 
QGP might be seen Qj. In the high-mass region at and 
above toj/^ the major effort has been the detection and 
understanding of J/\l/ suppression. So far, the experi- 
mental measurement of dilepton spectra at the CERN 
SPS has mainly been carried out by three collaborations: 
the CERES collaboration has specialized in dielectron 
spectra in the low-mass region the HELIOS-3 |t]] 

collaboration has measured dimuon spectra from thresh- 
old up to the J/fl> region, and the NA38/NA50 j§ col- 
laboration measures dimuon spectra in the intermediate- 
and high-mass regions. 

Recent observation of the enhancement of low-mass 
dileptons in central heavy-ion collisions over the proton- 
induced reactions by the CERES §,§ and the HELIOS-3 
JjJ collaborations has generated a great deal of theoreti- 
cal activity. The results from many groups with standard 
scenarios (i.e., using vacuum meson properties) are in re- 
markable agreement with each other, but in significant 
disagreement with the data: the experimental spectra 
in the mass region from 0.3-0.6 GeV are substantially 



underestimated ||. This has led to the suggestion of 
various medium effects that might be responsible for the 
observed enhancement PJTof. 

In the high-mass region around mj/^,, the J/^ sup- 
pression has been a subject of great interest, since it 
was first proposed as a signal of the deconfinement phase 
transition Various investigations show that up to 

central S+U collisions, the normal pre-resonance absorp- 
tion in nuclear matter is sufficient to account for the ob- 
served J/\l/ suppression. However, recent data from the 
NA50 collaboration for central Pb+Pb collisions show an 
additional strong 'anomalous' suppression which might 
indicate the onset of the color deconfinement . 

Another piece of interesting experimental data that has 
not received much theoretical attention is the dilepton 
spectrum in the intermediate-mass region between 1 and 
2.5 GeV. Both the HELIOS-3 and NA38/NA50 collab- 
orations have observed a significant enhancement of the 
dilepton yield in this mass region in central S-induced 
collisions, as compared to that in the proton-induced re- 
actions JjJH ■ Preliminary data from the NA50 collabora- 
tion also show significant enhancement in central Pb+Pb 
collisions (if (see also Ref. [[|). 

The intermediate-mass dilepton spectra in heavy-ion 
collisions are particularly useful for the search of the QGP 
Q . However, to extract from the measured dilepton spec- 
tra any information about the phase transition and the 
properties of the QGP, it is essential that the contribu- 
tions from the hadronic phase be precisely understood 
and carefully subtracted. It is the purpose of this work 
to calculate the dimuon spectra in central S+W colli- 
sions based on the relativistic transport model used in 
Ref. |Q for low -mass dilepton and photon production. 
The chief motivation for such a study is to understand the 
origin of the observed enhancement in the intermediate- 
mass region, and to see whether the enhancement can be 
explained by hadronic processes, or whether the forma- 
tion of the QGP needs to be invoked. 

Previous thermal rate calculations based on kinetic 
theory show that in the mass and temperature region 
relevant for this study, the following thermal processes 
(from the hadronic phase) are important: 7T7t — > 11, irp — > 
11, ttuj — > 11, nai — > 11, KK — > 11, and KK* + c.c — > 11 
]i~4] fl6| . Among them, the 7rai — > 11 has been found to 
be the most important, mainly because of its large cross 
section |l4|,[l6| (A similar conclusion has been drawn for 
thermal photon production jl7 18 1). In this work we shall 
verify quantitatively the previous estimates. 

To compare with experimental data, one needs a trans- 
port model that describes the dynamical evolution of the 
colliding system, and integrates the dilepton production 
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over the entire reaction volume and time. In heavy- 
ion collisions at CERN SPS energies, many hadrons 
are produced in the initial nucleon-nucleon interactions. 
This is usually modeled by the fragmentation of strings. 
One successful model for taking into account this non- 
equilibrium dynamics is the RQMD model fll9|| . As in 
Refs. |^,|l3|], we use as initial conditions the hadron abun- 
dance and distributions obtained from the string frag- 
mentation in RQMD. Further interactions and decays of 
these hadrons are then taken into account in a relativistic 
transport model. This model is found to provide a good 
description of hadronic observables in heavy-ion collisions 
at CERN SPS energies and does not rely on assumptions 
of thermal equilibrium . 

To calculate the dilepton spectra in heavy-ion colli- 
sions, we need to know the elementary cross sections for 
the secondary processes mentioned above. These pro- 
cesses can be classified into pseudoscalar-pseudoscalar 
(PP), vector-pseudoscalar (VP), vector- vector (VV), 
and axial-vector-pseudoscalar (AP) types. Relying on 
vector meson dominance (VMD), the first three types 
can be evaluated using effective hadronic Lagrangians 
0,^O|. Concerning the last case, there exist a number 
of models for npai dynamics. In Ref. pT| , a compara- 
tive study was carried out for both on-shell properties 
and dilepton production rates in those models. By us- 
ing the experimentally-constrained spectral function p2] ] 
it was found that the effective chiral Lagrangian of Ref. 
p3| in which the vector mesons are introduced as massive 
Yang-Mills fields provides satisfactory off-shell, as well as 
on-shell, properties for ixpa\ dynamics. The Lagrangian 
of Ref. | |23| will be used in this work in close conjunction 
with experimental data, as described below. 

The cross sections for these processes have a similar 
global structure. Here we take that for PP — > 11 as an 
example, 
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where k is the magnitude of the three-momentum of the 
pseudoscalar meson in the center-of-mass frame, M is 
the mass of the dilepton pair, mi the lepton mass, and 
\Fp(M)\ 2 is the electromagnetic form factor of the pro- 
cess. The cross sections for irp — > 11 and KK* +c.c. — > 11, 
which are of the VP type, have been studied in Ref. [ p^[ . 
The cross section for ttlo — * 11 has the same form as that 
for np — > 11, but with a different form factor. Finally the 
cross section for nai — > 11 is discussed below. 

We emphasize that these cross sections can be con- 
strained by a wealth of experimental data |2^] for the 
reverse process, e + e~ — > hadrons, via detailed-balance. 
One way this is achieved is by introducing vector meson 
dominance form factors which are fitted to the e + e~ — > 
hadrons data and then used in the meson- mesons 11 pro- 
cesses, as we will do in this work. The pion electromag- 
netic form factor is dominated by the p(770) meson, while 
that of the kaon is dominated by the 0(1020) meson. At 



large invariant masses, higher p-like resonances such as 
p(1450) were found to be important. High isoscalar vec- 
tor mesons such as w(1420) and 0(1680) play important 
roles in the electromagnetic form factors of irp — > 11 and 
KK* +c.c. — > 11 Jig] . Similarly, the electromagnetic form 
factor for ttoj — > 11 can be measured in e + e~ — > 7r°7r°7. 
The extraction of the form factor for ixa\ — > 11 from 
e + e~ — > 7T + 7r~7r + 7r - does involve some uncertainties con- 
cerning the p(1700) resonance. As this issue is not totally 
free of theoretical prejudice we will consider three sce- 
narios for the 7rai — > 11 form factor: one with and one 
without the p(1700) contribution, and one which uses a 
form factor determined from DM2 partial wave analy- 
sis data. The sensitivity of the final dimuon spectra to 
those cases will be assessed. A clear advantage of the 
e + e _ data is that its measurements cover the invariant 
mass region we are interested in here. Thus, no off-shell 
extrapolations are needed. Note that, however, the treat- 
ment of quantum interference remains a possible issue in 
our framework, see, e.g., Ref. |25|. 

For dilepton spectra with mass above 1 GeV, the con- 
tributions from charm meson decay and initial Drell-Yan 
processes begin to play a role. These hard processes, how- 
ever, scale almost linearly with the participant nucleon 
number, and can thus be extrapolated from the proton- 
proton and proton-nucleus collisions. Such a study has 
recently been carried out by Braun-Munzinger et al. [ p6| . 
The results for the central S+W collisions corresponding 
to the HELIOS-3 acceptance are shown in the left panel 
of Fig. 1 and are taken from Ref. ||. These, together 
with the dileptons from the decays of primary vector 
mesons, are collectively termed 'background'. It is seen 
that these background sources describe well the dimuon 
spectra in the p+W reactions, shown in the figure by 
solid circles. 

However, as can be seen from the figure, the sum 
of these background sources grossly underestimates the 
dimuon yield in central S+W collisions, shown in the fig- 
ure by open circles. Since the dimuon spectra are normal- 
ized by the measured charged particle multiplicity, this 
underestimation indicates additional sources of dilepton 
production in heavy-ion collisions. Some candidates are 
QGP and/or hadronic reactions. So the immediate next 
step is to check whether the contribution from the sec- 
ondary hadronic processes can explain this enhancement. 
For dilepton spectra at low invariant masses, it is well 
known that the nir annihilation plays an extremely im- 
portant role in heavy- ion collisions. It is also expected 
that the other secondary processes will play a role in the 
dilepton spectra in the intermediate mass region. 

The contributions from the thermal processes outlined 
above are shown in the middle panel of Fig. 1. These 
are obtained in the relativistic transport model of Refs. 
HlH, including the HELIOS-3 acceptances, mass reso- 
lution, and normalization [Q. It is seen that the nai 
process is by far the most important source for dimuon 
yields in this mass region, as was found in thermal rate 
estimates. The nui process also play some role in the 
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entire intermediate mass region, while the contributions 
from 7T7T, 7T/9 and KK are important around 1 GeV. We 
have also verified that the AV (where A = ai) contribu- 
tions are suppressed P7J. 
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FIG. 1. Left panel: comparison of background and pre- 
liminary data in p+W and S+W collisions. Middle panel: 
contributions of various secondary processes of dimuon pro- 
duction in central S+W collisions. Right panel: comparison 
of the sum of the background and secondary contributions 
with the preliminary data in central S+W collisions. 



In the right panel of Fig. 1, we add the secondary 
contributions to the background, and compare again 
with the HELIOS-3 data for central S+W collisions. 
It is seen that the data can be adequately reproduced. 
This highlights for the first time the importance of the 
secondary processes for the intermediate-mass dilepton 
spectra in heavy-ion collisions even though the evalua- 
tion of the background itself is not entirely free of com- 
plications. This is an important step forward in the use 
of intermediate-mass dilepton spectra as a probe of the 
phase transition and QGP formation. Although the cur- 
rent data do not show any necessity to invoke the QGP 
formation in S-induced reactions, consistent with some 
conclusions from J /"J physics, we believe that the obser- 
vation that the secondary processes do play a significant 
role in the intermediate-mass dilepton spectra is inter- 
esting and important. We note in passing that the slight 
change of the slope observed in the experimental data 
corresponds in our interpretation to a crossover between 
the secondary processes and the (Drell-Yan) background. 
This signals a passage from hard to soft physics and could 
lead to interesting developments. 

In the previous calculation we assumed that all the 



ceeds through the nai intermediate state. We also did 
a calculation in which the irai form factor contains only 
the normal p(770). The results are shown in the left 
panel of Fig. 2 by the dotted curve. Finally, to complete 
our survey of possible constraints and uncertainties in 
■Ka\ — > e + e~ cross sections, we did a calculation in which 
this cross section is obtained from the e + e~ — > nai cross 
section determined by the DM2 collaboration in partial 
wave analysis (PWA) pq |. The results are shown by 
dashed curve in Fig. 2. The region between the solid 
and dashed curves thus reflect the uncertainty for heavy- 
ion collisions due to our limited knowledge of the nai 
cross section. This area is not unreasonably large. From 
a formal point of view, it is fair to say that no strong 
evidence currently exists coupling the p(1700) to a irai 
state ]2!| , even though better 4-7T data could help resolve 
this issue along with others of interference and nhi con- 
tribution 
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FIG. 2. Left panel: comparison of dimuon spectra obtained 
with and without the p(1700) in the ivai form factor, with 
those obtained with the DM2 data. Right panel: compari- 
son of dimuon spectra with im-medium and vacuum meson 
masses. 



Another topic to be adressed here is the effect of drop- 
ping vector meson masses on the entire dimuon spec- 
tra from threshold to about 2.5 GeV. In Ref. |J it was 
shown that the enhancement of low-mass dileptons could 
be interpreted as a signature of vector meson masses de- 
creasing with increasing density and temperature. This 
should affect the dilepton spectra in the intermediate- 
mass region, mainly through two effects. One is the 
change of the invariant energy spectra of these secondary 
meson pairs. The second effect enters through the mod- 
ification of the electromagnetic form factor. Since we 
can only conjecture how the masses of the higher vec- 
tor resonances change with density and temperature, 
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we shall assume for simplicity that they experience the 
same amount of scalar field as the "common" rho me- 
son, namely, rriy V , = myy — 2/3g a (a) The results 
of this calculation are shown in the right panel of Fig. 
2. Below 1.1 GeV and especially from 0.4 to 0.6 GeV, 
the agreement with the experimental data is much bet- 
ter when the dropping vector meson mass scenario is in- 
troduced, as was already shown in Ref. |J. At higher 
invariant masses the dropping mass scenarios somewhat 
underestimates the experimental data, admittedly not by 
a large amount. For completeness, however, we have to 
state that there might be additional contributions from, 
e.g., secondary Drell-Yan processes |30j that were not in- 
cluded in this study. Furthermore, as we progress higher 
in invariant mass, the role of baryons has to be carefully 
assessed. So far, the baryons seem to play little role in 
the overall dilepton yield J3l| . This statement was made 
p2"| for masses below 1 GeV, and is being extended to 
the intermediate-mass region in f33f| . See however Ref. 

for a contrasting viewpoint on the role of baryons. 
Finally, the dropping mass curve was obtained using the 
topmost curve of the left panel as a starting point. 

So far, collision broadening is not explicitly included 
in our calculation. However, since the "common" rho 
meson is treated as a dynamical particle in our trans- 
port model, its collisions with mesons and baryons are 
included. This shall partially reflect its collision broad- 
ening in hot and dense matter fsif . The higher reso- 
nances, such as p(1450) and p(1700), are not treated dy- 
namically. Their effects on dilepton production are in- 
cluded through electromagnetic form factors. Therefore, 
collision broadening is not considered explicitly for these 
higher resonances. However, Since they typically have a 
natural width of 200-300 MeV, a moderate broadening 
width shall not affect our results dramatically. This will 
be dealt with quantitatively in a more elaborate way in 
the future. 

In summary, we have analysed the recent HELIOS-3 
data on intermediate-mass dilepton production in heavy- 
ion collisions at CERN SPS energies using a relativistic 
transport model, with the elementary dilepton produc- 
tion cross sections constrained by the e + e~ annihilation 
data. We have shown the importance of secondary pro- 
cesses for the dilepton production in heavy-ion collisions 
in this mass region. 

The current investigation can be extended to higher 
incident energies, such as those of the RHIC collider, 
by combining the cross sections (or thermal rates) ob- 
tained in this study with, e.g., hydrodynamical models 
for the evolution of heavy-ion collisions at the RHIC en- 
ergies. This kind of study is useful for the determination 
of hadronic sources in the dilepton spectra, and for the 
clear identification of the dilepton yield from the QGP. 
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